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Abstract 

Among some 400 million KsKl pairs produced in e~^e~ annihilations at DA<I>NE, 
~6500 each of i^^— >7r"'"e~i> and K5— >7r~e+i^ decays have been observed with the 
KLOE detector. From these, the ratio r{Ks 'Kev)/V{Ks tt+tt") = (10.19 ± 
0.13) X 10~^ is obtained, improving the accuracy on BR(i^5 nev) by a factor 
of four and providing the most precise test of the A5 = AQ rule. From the partial 
width T{Ks — > vrez^), a value for f^ (0) x V^s is obtained that is in agreement 
with unitarity of the quark-mixing matrix. The lepton charge asymmetry As = 
(1.5±9.6stat±2.9syst) X 10~^ is compatible with the requirements of CPT invariance. 
The form-factor slope agrees with recent results from semileptonic Kl and 
decays. These are the first measurements of the charge asymmetry and form- factor 
slope for semileptonic Kg decays. 



1 Introduction 

Semileptonic kaon decays provide at present the best way to learn about s, u 
quark couplings and allow tests of many fundamental aspects of the Standard 
Model (SM). Only the vector part of the weak current has a non- vanishing 
matrix element between a kaon and a pion. The vector current is "almost" 
conserved. For a vector interaction, there are no S'f/(3)-breaking corrections 
to first order in the s-d mass difference [1], making calculations of hadronic 
matrix elements more reliable. Therefore, the CKM matrix element V^g can be 
accurately extracted from the measurement of the semileptonic decay widths 
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and the most precise test of unitarity of the CKM matrix can be obtained 
from the first-row constraint: 1 — A~ |KidP + |KtsP- Using V^d from 0"^ — ^ O"*" 
nuclear beta decays, a test of the expectation A = with a precision of one 
part per mil can be performed. 

At a factory very large samples of tagged, monochromatic Ks mesons are 
available. We have isolated a very pure sample of ~13 000 semileptonic Ks 
decay events and measured for the first time the partial decay rates for transi- 
tions to final states of each charge, T{Ks^e'^7r~i') and r{Ks-^e~7r~^iy), and 
the charge asymmetry 

The comparison of with the asymmetry Al for Kl decays allows tests of the 
CP and CPT symmetries. Comparison of the Kg and widths T{Ks 
Treu) and T{Kl Tieu) allows a test of the validity of the AS* = AQ rule. 
Assuming CPT invariance, As = A^ = 2Ree ~ 3 x 10"^, where e gives the 
CP impurity of the Ks, Kl mass eigenstates due to CP violation in K ^ K 
AS — 2 transitions. The difference between the charge asymmetries, 

As-AL^4:{ReS + Rex^), (2) 

signals CPT violation either in the mass matrix {6 term), or in the decay 
amplitudes with ASy^AQ {x- term). The sum of the asymmetries, 

As + AL^4.(Ree-Rey), (3) 

is related to CP violation in the mass matrix (e term) and to CPT violation in 
the AS = AQ decay amplitude {y term). Finally, the validity of the AS = AQ 
rule in CPT-conserving transitions can be tested through the quantity: 

_ 1 r{Ks ^ TTCiy) - r{KL ^ Treiy) 
+ 2 r{Ks ^ neu) + r{KL ^ neu) ' ^ ' 



Writing the K^ and K^ decay amplitudes for final states of each charge as 
A± — A{K^ ^ e^7r^i'{iy)) and A± — A{K°^e^7r^i'{iy)), the above parameters 
are defined as follows: 




(5a) 



At + A. 
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ImMi2 - iImri2/2 



e 



ms -rriL- i{l/rs - 1/tl)/2' 

1 Mn - M22 - ^(^n - r22)/2 

2 ms -rriL- K'^/ts - 1/tl)/2' 



(5b) 
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where Mij and Fjj are the elements of the mass and decay matrices describing 
the time evolution of the neutral kaon system, and ms,L and ts,l are respec- 
tively the masses and lifetimes for Ks,l- 

The value ol Ai^ is known at present with an accuracy of 10~^ [2], while 
has never yet been measured. At present, the most precise test of CPT con- 
servation comes from the CPLEAR experiment [3]: they find Re 5 and Rea;_ 
to be compatible with zero, with accuracies of 3x10^'^ and 10^^, respectively. 
The present value of Re|/, obtained from unitarity [4], is compatible with zero 
to within 3 x 10~^. 

In the SM, Rea;+ is on the order of Gpfn^^ 10^^, being due to second-order 
weak transitions. At present, the most precise test of the A5' = AQ rule comes 
from an analysis of the time distribution of strangeness-tagged semileptonic 
kaon decays at CPLEAR [5] . They found Re x+ to be compatible with zero to 
within 6x10-^. 

The most precise previous measurement of ^Ii{Ks — > Tiev) was obtained by 
KLOE using ~ 2C^b~^ of data collected in 2000 and has a fractional accuracy 
of 5.4% [6]. The present result is based on the analysis of 410 pb~^ of data 
and improves on the total error by a factor of four, to 1.3%. 



2 Measurement method 

We measure Ks branching ratios using kaons from KsK^ decays. The 
data were collected with the KLOE detector at DA$NE, the Frascati 0- 
factory. DA$NE is an e^e~ collider that operates at a center of mass energy 
of ~1020MeV, the mass of the meson. Equal-energy positron and electron 
beams collide at an angle of tt — 25 mrad, producing mesons with a small 
momentum in the horizontal plane: 13MeV. mesons decay ~34% of the 
time into neutral kaons. K^s and Ks^s have mean decay paths of AL~350cm 
and A5~0.6cm, respectively. 

The KLOE detector consists of a large cylindrical drift chamber surrounded by 
a lead/scintillating- fiber sampling calorimeter. A superconducting coil outside 
the calorimeter provides a 0.52 T field. The drift chamber [7], which is 4m in 
diameter and 3.3 m long, has 12 582 cells strung in all-stereo geometry. The 
chamber shell is made of a carbon-fiber/epoxy composite. The chamber is filled 
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with a 90% He, 10% iC4Hio mixture. These features maximize transparency to 
photons and reduce K^^Ks regeneration. The spatial resolutions are a^y ~ 
150 /xm and cr^^ 2 mm. The momentum resolution is (Jp^)/p±<OA%. Vertices 
are reconstructed with a spatial resolution of ~ 3 mm. The calorimeter [8] 
is divided into a barrel and two endcaps and covers 98% of the solid angle. 
The energy resolution is ue/E — 5.7%/ ^E{GeV) and the timing resolution 

is (7t = 57ps/y£^(GeV) © 100 ps [9]. The trigger used for the present analysis 
relies entirely on calorimeter information [10]. Two local energy deposits above 
threshold (50 MeV on the barrel, 150 MeV on the endcaps) are required. The 
trigger time has a large spread with respect to the bunch-crossing period. 
However, it is synchronized with the machine RF divided by 4, 

-^sync ^ 10.8 ns, 

with an accuracy of 50 ps. As a result, the time Tq of the bunch crossing 
producing an event, which is determined after event reconstruction, is known 
up to an integer multiple of the bunch-crossing time, 7bunch~ 2.7ns. 

The main advantage of studying kaons at a factory is that K^s and Ks^s 
are produced nearly back-to-back in the laboratory so that detection of a Kl 
meson tags the production of a Ks meson and gives its direction and mo- 
mentum. The contamination from KiKcj and KsK^ final states is negligible 
for our measurement [11, 12]. Since the branching ratio for Ks 7r~^7r~ is 
known with an accuracy of ~ 0.1% [13,14], the Ks — > Treiy branching ratio 
is evaluated by normalizing the number of signal events, separately for each 
charge state, to the number of Ks tt'^tt' events in the same data set. This 
allows cancellation of the uncertainties arising from the integrated luminosity, 
the 4> production cross section, and the tagging efficiency. The measurement is 
based on an integrated luminosity of 410 pb~^ at the peak collected during 
two distinct data-taking periods in the years 2001 and 2002, corresponding to 
~1.3 X lO'^ produced 0-mesons. Since the machine conditions were different 
during the two periods, we have measured the branching ratios separately for 
each data set. Our final results are based on the averages of these measure- 
ments. 



3 Selection criteria 



About half of the Kl mesons reach the calorimeter, where most interact. Such 
an interaction is called a K^ crash in the following. A K^ crash is identified 
as a local energy deposit with E > 200 MeV and a time of flight corresponding 
to a low velocity: /3~ 0.216. The coordinates of the energy deposit determine 
the Kl direction to within ~20mrad, as well as the momentum Pl, which is 
weakly dependent on the Kl direction because of the motion of the meson. 
A Kl crash thus tags the production of a Ks of momentum Ps = p</> — Pl- 
Ks mesons are tagged with an overall efficiency of ~ 20%. Both Ks — > neu 
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Year 2001 Year 2002 

Luminosity, pb~^ 152 ± 1 286 ± 2 



Channel Number of selected events 



Ks-^ 


■ TT+TT" 


13 056 500 


22 840 700 


Ks^ 




2387 ± 52stat i 22syst 


4238 ± 69stat ± 55syst 


Ks^ 




2541 ± 52stat ± 24syst 


4446 ± 71stat ± 69syst 



Table 1 

Number of selected Ks tt+tt" and Ks ireu decays for the 2001 and 2002 data 
sets. 

and Ks — * tt+tt" decays are selected from this tagged sample. Event selection 
consists of fiducial cuts, particle identification by time of flight, and kinematic 
closure. 

Identification of a Ks — > tt+tt" decay requires two tracks of opposite curvature. 
The tracks must extrapolate to the interaction point (IP) to within a few 
centimeters. The reconstructed momenta and polar angles must lie in the 
intervals 120 MeV <p < 300 MeV and 30° < ^ < 150°. A cut in {p±,p\\) selects 
non-spiralling tracks. The numbers of Ks tt+tt" events found in each data 
set are shown in Tab. 1. Contamination due to Ks decays other than Ks — > 
7r"'"7r~ is at the per-mil level and is estimated from Monte Carlo (MC). 

Identification of a Ks — > vrez/ event also begins with the requirement of two 

tracks of opposite curvature. The tracks must extrapolate and form a vertex 
close to the IP. The invariant mass Af^^ of the pair calculated assuming both 
tracks are pions must be smaller than 490 MeV. This rejects ~ 95% of the 
7r+7r~ decays and ~ 10% of the signal events. 

Wc discriminate between electron and pion tracks by time of fiight (TOF). 
The tracks are therefore required to be associated with calorimeter energy 
clusters. For each track, we compute the difference St{m) = ta — L/c/3{m) 
using the cluster time ^ci and the track length L. The velocity P is computed 
from the track momentum for each mass hypothesis, m = me and m = mT^. In 
order to avoid uncertainties due to the determination of To (the time of the 
bunch crossing producing the event), we make cuts on the two-track difference 

where the mass hypothesis ma(b) is used for track 1(2). This difference is zero 
for the correct mass assignments. First, Ks — > tt+tt^ events are rejected by 
requiring |d(5t,^7r| > 1.7 ns. Then, the differences d^j^^e and dSt^^Tr are calculated 
for surviving events. The scatter plot of the two variables is shown in Fig. 1 for 
Monte Carlo events. The cuts applied on these time differences for the selection 
of Ks — > TTCU events are illustrated in the figure: |d(5t^,re| < 1.4ns, d(5(,e7r >3.2 ns; 
or |d5f,e7r| < 1.4 ns, d5t,7re> 3.2 ns. After these TOF requirements, particle types 
and charges for signal events can be assigned very precisely: the probability of 
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Fig. 1. Scatter plot of the time differences ddt^ne vs ddt^en for Monte Carlo events, 
for all Ks decays (left) and for Ks — > nev decays (right). 

misidentifying a ir^e^u event as 7r~e"^z/ or vice versa is negligible. These cuts 
reject ~ 90% of the background events, while the efficiency for the signal is 
~85%. 



Once particle identification has been performed, we reevaluate the time differ- 
ences 5t{'m), this time using for each track the mass assignment known from 
the cut on dSt^eir and subtracting the Tq of the event. For the Tq determina- 
tion, the bunch crossing producing the event is evaluated as the integer part 
of the ratio [5t{e) + 5t(7r)]/2Tbunch- We apply another TOF cut by selecting 
the events within the circle in the 5t{e)-5t{T^) plane, as shown in Fig. 2 for MC 
events. This cut improves the background rejection by a factor of five, while 
eliminating 8% of the signal events at this analysis stage. 



-2 - 




-2 








-2 



-1 







Fig. 2. Scatter plot of the time differences (5j(7r) vs 5t{e) for tt and e mass assignments 
for Monte Carlo events, for all Ks decays (left) and for Ks Trez/ decays (right). 
Events within the circles are retained. 
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A powerful discriminating variable is the difference between the missing en- 
ergy and momentum, AE^^e = -Emiss ~ Pmissj which is evaluated using the Ks 
momentum known from the direction. For ttci' decays, £^miss and Pmiss are 
the neutrino energy and momentum, and are equal. The distribution of AE^^. 
is shown in Fig. 3 after TOF cuts are imposed for 7r~e'^u (left panel) and for 
7r'^e~P (right panel) candidate events. A clear peak around zero is evident and 
corresponds to a clean signal for Ks Trei/. 

The residual background is dominated by Ks — > 7r''"7r~(7) decays. Events 
with AEj^g > lOMeV are mostly due to cases in which one pion decays to a 
muon before entering the tracking volume ("vr/i" events), in which the track 
identified as electron by TOF is badly reconstructed ("TTbadTr" events), or in 
which the radiated photon has an energy in the Ks frame above 7MeV, thus 
shifting Mj^T^ below 490 MeV and £'miss toward positive values ( "7r7r7" events) . 
Events with AE^^^ < —10 MeV are mostly "tt//" or "TTbadTr" events, or are due to 
cases in which the track identified as the pion by TOF is badly reconstructed 
("TTTTbad" events). 

We discriminate between signal and residual background events by means of 
5 kinematic variables: AE'^g itself; the difference dpcA between the impact 
parameters of the two tracks with the IP; the difference AE'^tt in the tt-tt 
hypothesis; the squared mass M^j.y^{e) of track 1(2) when it is identified as 
electron from TOF, calculated assuming that ps — P2{i) is the momentum of 
an undetected pion and that {ps — Pi — ^2)^ = 0; the energy E*^^^-^ of the track 
identified as a tt (e) from TOF, calculated in the Ks rest frame using the pion 
mass hypothesis. 

Except for "7r7r7" events, all of the background categories arc characterized 
by poor vertex reconstruction quality, leading to a broad distribution of dpcA 
as shown in Fig. 4, top left. In contrast, signal and "71717" events are peaked 
around zero. We discriminate "tt/x" events by the M^^^{e) variable, which peaks 
around m^, and "71717" events by the AEj^j^ variable, which peaks around zero. 
Finally, well reconstructed pion tracks from Ks 7r~^7i~ decays arc identified 
by the value of E*, which peaks around niK/^, allowing us to recognize 7rbad'?'" 
or '^'^bad events. 

The number of events due to signal and to each of the background categories 
are evaluated through a global binncd-likelihood fit using the above variables. 
Each event is assigned to one of five regions in the AE'^g-dpcA plane illustrated 
in Fig. 5. For each region, we use for the fit one of the variables defined above. 
The choice of the regions and the assignment of the fit variables to each region 
ensure good separation between each component in turn and all the others. 
In each region, the data are fit with the sum of the MC distributions in the 
appropriate variable for signal events and for events from each background 
source. The free parameters are the signal and background normalizations. 
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Fig. 3. E'miss — f'miss spectrum for events selected as vr^e^z/ (left panel) and 7r~^e~i' 
(right panel). Filled dots represent data from the entire data set; the solid line is 
the result of a fit varying the normalization of MC distributions for signal (light 
gray) and background (dark gray), which are also shown. 



For each source, the same normalization parameter is used in all fit regions. 



The result of the fit is shown as the solid hne in the distributions of Figs. 3 

and 4. The MC simulations of Ks vrez/ and Ks ti^ti~ decays include 
photon radiation in the final state [15]. If the effect of radiation were not 
taken into account, the result for the branching ratio would decrease by a few 
percent. 

We perform two independent fits, one for each charge state. The estimated 
numbers of signal events are shown in Tab. 1. The quoted statistical errors 
include the contributions from fluctuations in the signal statistics, from the 
background subtraction, and from the finite MC statistics [16]. 

The systematic errors include the contribution from uncertainty in the shape of 
the signal distributions. In particular, we have studied in detail the reliability 
of the MC in reproducing the distribution of AE'^e- We have compared data 
and MC resolutions obtained from samples of Kg — > tt+tt" events tagged by 
the Kl crash, both for the momentum of each track and for the AE^^j^ variable. 
From these studies, we have extracted corrections to the MC resolutions for 
the tracking momentum and the polar and azimuthal angles for the Kl crash. 
After these corrections are applied, good agreement is observed for the core 
of the AE'^e distribution for signal events. In order to study the rehabihty 
of the simulation for the tails, an additional variable xpiD has been used to 
discriminate signal events, xpio identifies e"^ tracks on the basis of the spatial 
distribution of the energy deposition in the EmC, and is independent from 
the track momenta used to obtain the fit variables. An alternative estimate 
of the number of signal events in each fit region is obtained from the xpiD 



9 




-10 10 20 
Mttk(e)/10^ (MeV^) 
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Fig. 4. Distributions for Tr'^e~i/ candidate events, for data (dots) and Monte Carlo 
events (solid line). Top left: dpcA variable, peaking around zero for vertices of 
good quality (signal and inr'y events); top right: AEj^t^, peaking around zero for 
7r7r7 events; bottom left: M^^-^^{e), peaking around for tt/j, events; bottom right: 
E*, peaking around 771^-/2 for TTTTbad events. In each plot, signal (light gray) and 
background (dark gray) contributions are also shown. 

distribution, which is reliably reproduced by the MC simulation. Using this 
method, a significant difference in the number of signal events is observed only 
in regions 1 and 2, from which a shape correction is defined. This corresponds 
to a 0.5% correction to the final result, which is also taken as the corresponding 
systematic error. 



The cut on the minimum cluster energy required for KL-crash tagging dramat- 
ically affects the fraction of signal events in the tails of the AE'^e distribution: 

the looser the cut, the worse the resolution on the Ks momentum evaluated 
from the direction. In order to check the robustness of the signal extrac- 
tion, we have compared results obtained using three values for the minimum 
Kl cluster energy: 125 MeV, 200 MeV, and 300 MeV. The results obtained us- 
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Fig. 5. Definition of the fit regions in the Ai?7re-dpcA plane. In each region, the 
variable used for the fit is also specified. 

ing independent samples arc compatible with each other. The intermediate 
cut gives the minimum total uncertainty. 

A detailed description of the fit procedure is given in Ref. 9. 



4 Efficiency estimates 



For both Ks tt+tt^ (normalization) and Ks vrez/ (signal) events, we 
estimate the corrections for the geometrical acceptance and tagging and selec- 
tion inefficiencies with the MC simulation. In order to account for data-MC 
differences, we weight each event with the ratio of the data and MC tracking 
efficiencies extracted from various control samples. We evaluate the probabil- 
ities for finding EmC clusters from Ks daughter particles and for satisfying 
the trigger conditions by combining data-extracted efficiencies, parametrized 
in terms of track momenta, with MC kinematics. For Ks — > Treu events, we 
also use prompt Kl — > Treu decays tagged by Ks — > tt'^tt' as a data control 
sample for efficiency evaluation. With this method, we obtain alternative es- 
timates for the trigger and cluster efficiencies and evaluate the corrections for 
vertex reconstruction and vr-e TOF identification inefficiencies. These methods 
are described in detail in Refs. 9 and 17. 

For Ks vrez/ decays, the efficiencies are determined separately for each 
charge state and for the 2001 and 2002 data. We summarize the results for 
the overall efficiencies, given the tag requirement, in Tab. 2. The differences 
between the efficiencies for the two charge states arise from the different re- 
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Kg decay 


Selection efficiency 
Year 2001 Year 2002 




0.5954 ± 0.0004stat ± O.OOlOgyst 0.6035 ± 0.0004stat ± O.OOlOgygt 




0.2139 ± 0.0019stat ± O.OOUgygt 0.2197 ± 0.0012gtat ± 0.0021syst 




0.2252 ± O.OOlGgtat ± O.OOODgygt 0.2328 ± O.OOllgtat ± O.OOllgyg^ 



Table 2 

Selection efficiencies for Kg tt~^tt~ and Kg 'Keu decays, for the 2001 and 2002 
data sets, given the i^L-crash tag. 



sponse of the calorimeter to tt"*" and 7r~, influencing the cluster, trigger, and 
TOF efficiencies. The uncertainties on the tracking and cluster/trigger efficien- 
cies contribute approximately equally to the systematic errors on the overall 
efficiencies. A variation in the cluster/trigger efficiencies between 2001 and 
2002 is reflected in the values for the overall efficiencies. The corresponding 
systematic errors have been estimated from the comparison of the results ob- 
tained using prompt Tieu decays and using single-particle weights to 
correct the MC. The difference between the results obtained with the two 
methods is larger for 7r~e+i/ events, and for the 2002 data. 

In principle, the i^^L-crash identification efficiency cancels out in the ratio of 
the number of selected Kg — > neu and Kg — > tt+tt^ events. In practice, since 
the event Tq is obtained from the Kg and the Kl is recognized by its time 
of flight, there is a small dependence of the fr^-crash identification efficiency 
on the Kg decay mode. A correction for this effect is obtained by studying 
the accuracy of the To determination in each case [9, 17, 18]. The ratio i?tag of 
the tagging efficiencies for Kg Treu and Kg —>■ tt+tt" is found to differ from 
unity by ~ 1%. This effect is included in the efficiency values shown in Tab. 2. 



5 Results 



For each charge state and for the data set for each year, we obtain the ratio of 
BR's by normalizing the number of signal events to the number of Kg — > tt+tt" 
events and correcting for the overall selection efficiencies: 

r{Kg^Tr+7r-) Ninn) e±t' 

The results for the BR's from the data sets for each year are compatible with 
probabilities greater than 50%. Averaging the results obtained for each data 
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set, we obtain the following results: 



(5.099 ± 0.082,tat ± 0.039syst) x 10"^ 

(5.083 ± 0.073stat ± 0.042,y,t) x 10"" (6) 

(10.19 ± O.llstat ± O.OTsyst) X 10-^ 

To obtain the value of the ratio BR{Ks — > 7rei')/BR{Ks — > tt+tt") we take 
into account the correlation between the values measured for the two charge 
modes. This correlation arises from uncertainties on the shapes of the signal 
distributions in the fit variables that are common to both charge states; the 
correlation parameter is 13%. 

The charge asymmetry of Eq. (1) is given by: 

, ^ 7V(7r-e+i/)/e+t - N {tt^ e- u) / e^,^ 
^ N{n-e+u)/e+t + N{n+e-u)/e^,, ' 

Combining the results for all data, we obtain: 

^5 = (1.5 ± 9.6stat ± 2.9syst) X 10-^ 

In order to perform a stability test, we have divided the entire data set into 17 
samples and performed the analysis individually for each sample. Values of 
corresponding to probabilities above 50% are observed for all of the measured 
quantities [9]. 

The various contributions to the total fractional error on BR{Ks Treu) and 
to the total error on As are listed in Tabs. 3 and 4. For the measurement of the 
BR, the uncertainty on the signal count from fit systematics is the dominant 
contribution to the total systematic error. 

For the purposes of measuring the dependence of the form factor /+(t)//+(0) 
on the 4- momentum transfer squared t — {Pk — P-kT, the elimination of back- 
ground from the sample while preserving statistics is a more important consid- 
eration than the understanding of the selection efficiencies at the level required 
in the branching-ratio analysis. We therefore use slightly different selection cri- 
teria to isolate a clean sample of Ks irei' events: we require | Aii^Trel < 10 MeV 
and cut on the quality of the Ks vertex. In order to limit loss of statistics, we 
loosen the energy requirement on the crash to 125 MeV. We select about 
15 000 signal events, combining the data for the two years and charge modes. 
With this selection, the bacgkround contamination is reduced to 0.7%. Be- 
cause of the limited statistics, we only measure the slope parameter of the 
form factor in the hnear approximation, /+(t)//+(0) = 1 + A+t/m^. More 



_ T{Ks — > TT e'^u) 

" r{Ks TT+TT-) 

, _ r(A-,, ^ n+e-u) 

r{Ks^ TT+TT-) 

^ ^ T{Ks TTeu) 
^ r(Ks 7r+7r~) 
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Fractional 


error (10 ^) 




Statistical 


Systematic 


Statistics of Kg —>■ T^ev 


9.1 




Statistics of Kg ■n^ir" 


0.1 




Preselection efficiency, Kg ttcu 


1.5 


2.9 


Trigger efficiency, Kg ireu 


0.2 


0.3 


TOF efficiency, Kg —>■ Trev 


2.3 




Fit systematics, Kg irev 




6.2 


Preselection efficiency. Kg —>■ tt^tt^ 


0.3 


1.6 


Trigger efficiency, Kg — tt+tt" 


0.1 


0.8 


Ratio of tagging efficiencies 


5.0 




Ratio of cosmic veto inefficiencies 


1.0 




Total 


10.8 


7.1 



Total fractional error 12.9 



Table 3 

Contributions to the fractional error on BR{Kg —>■ ■Kei')/BR{Kg — > 7r"'"7r~). 



Error (lO"-'') 
Statistical Systematic 



Statistics of Kg irev 


9.1 




Preselection efficiency. Kg — > nev 


1.5 


2.9 


Trigger efficiency. Kg — > Treu 


0.1 


0.3 


TOF efficiency, Kg — >■ weu 


2.3 




Fit systematics, Kg — >■ weu 




0.4 


Tagging efficiencies 


0.4 




Cosmic veto inefficiencies 


1.0 




Total 


9.6 


2.9 



Total error 



Table 4 

Contributions to the absolute error on Ag. 



10.0 



precisely, we fit the ratio of data and MC distributions in t/m 
function: 

I I —I— yi I / / ■/ 1^ ^ 

F{t) = A X 



with the 



1 + X+t/ml, 



1 + A_ 



-,MC 



t/ml^ 



where A and A+ are the free parameters of the fit, and X+^mc = 0.03 is the 
value of the slope used in the MC generation. Effects from the finite resolution 
on t are negligible with respect to the statistical error and are ignored. We find 
A+ = (33.9±4.1) X 10^^ with x^/dof = 12.9/11, corresponding to a probability 
P(x^) — 30%. This result is in reasonable agreement with the value of A+ for 
semileptonic Kl and decays, (28.82 ± 0.34) x 10~^, from the average of 
results from KTeV [19], ISTRA+ [20], NA48 [21], and KLOE. 
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6 Interpretation of the results 



6.1 Determination of absolute BR 's 



In order to evaluate the BR's for the semileptonic modes, we combine the ratios 
of BR's measured for each charge [Eq. (6)] with the most precise measurement 
of the ratio 

T(Ks 



R 



TT ' TT 



2.236 ±0.015, 



(7) 



r{Ks TT^TrO) 

which was also obtained at KLOE [23]. The only remaining mode with a 

BR large enough to measurably affect the constraint J2 / BR(i^5 ^ /) = 1 
is K^s] the BR's for all other channels sum up to ~ 10~^. Assuming lepton 
universality, we have 



(8) 



where are mode-dependent long-distance radiative corrections and /^'^ are 
decay phase-space integrals. Using Ik/^k — 0.6622(18) from KTeV [24] and 
(1 + + ^k) = 1.0058(10) from Ref. 25, we get r^e = 0.6660(19). We 

evaluate the four main BR's of the Ks from 



BR{Ks ^ i) = 



T{Ks i)/TiKs TT+TT-) 

l + l/i?+(i?e+ + i?e-)(l + V)' 



(9) 



where i — TT'^ir , 7r°7r°, tt e'^u, ir'^e v. We find: 



BR{Ks ttV) 
BR(ii'5 ^ TT-e+i/) 



(69.02 ±0.14) X 10-2 
(30.87 ±0.14) X 10-2 
(3.519 ± 0.063) X 10-^ 
(3.508 ±0.058) X 10"^ 



(10) 



The correlation matrix {SiSj)/J{5f){Sj) is: 



TT ' e u 



TT+TT 



tt'^ttO 



TT e+i^ TT"'"e u 



1 -0.9999 0.1106 0.1196 ^ 
-0.9999 1 -0.1187 -0.1272 
0.1106 -0.1187 1 0.1445 
0.1196 -0.1272 0.1445 1 y 



:iii 
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The contribution from the error on r^e is included in the systematic errors. 
Taking correlations into account, we have: 

BR{Ks neu) = (7.028 ± 0.092) x 10"^ (12) 

6.2 Test of the AS = AQ rule 

From the total BR we test the validity of the A^* = AQ rule in CPT- 
conserving transitions [Eq. (4)]. We use the following values for the Ks and 
Kl lifetimes: ts = 0.08958(6) ns from the PDG [13] and tl = 50.84(23) ns 
from recent measurements from KLOE [26,27]. For BK(Kl — > neu), we use 
the value from KLOE, 0.4007(15) [26]. We obtain: 

Rex+ = (-1.2±3.6) X 10-^ (13) 

The error on this value represents an improvement by almost a factor of 
two with respect to the most precise previous measurement, that from the 
CPLEAR experiment [5]. 

6.3 Test of the CPT symmetry 

From the sum and difference of the Kl and Ks charge asymmetries one can 
test for possible violations of the CPT symmetry, either in the decay ampli- 
tudes or in the mass matrix [Eqs. (2) and (3)]. Using Al — (3.34 ± 0.07) x 
10-^ [13], we obtain from Eq. (2) 

Re a;_ + Re 5 = (-0.5 ± 2.5) x 10"^ (14) 

Current knowledge of these two parameters is dominated by results from 
CPLEAR [3]: the error on Re 5 is 3 x 10"*^ and that on Rex- is 10"^. Us- 
ing Re5 = (3.0 ± 3.3stat ± 0.6syst) x 10"^ from CPLEAR, we obtain: 

Rex_ = (-0.8 ± 2.5) x 10-^ (15) 

thus improving on the error of Re x_ by a factor of five. 

From Eq. (3) we obtain 

Ree-Rey = (1.2±2.5) X 10-^ (16) 

We calculate Re e using values from Ref. 13 that are obtained without assum- 
ing CPT conservation: Ree = |e| x cos^g = (1.62 ± 0.04) x 10~^. Subtracting 
this value from Eq. (16), we find 

Rey = (0.4±2.5) X 10-^ (17) 
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which has precision comparable to that (3x10 ^) obtained from the unitarity 
relation by CPLEAR [4]. 



6.4 Determination ofV^ 



The value of Vus can be extracted from the measurement of BR(A'5 nev) 
and from the Ks lifetime, r^: 



tsGIMI,S,^Ik{X+,0) 



1 + 



(18) 



where f^°^ (0) is the vector form factor at zero momentum transfer and 
Ik{X+, 0) is the result of the phase space integration after factoring out f^"'^ (0); 
both quantities are evaluated in absence of radiative corrections. The radia- 
tive corrections [25, 28] for the form factor and the phase-space integral are 
included via the parameter 5^ = (0.55 ± 0.10) xlO^^ [28]. The short-distance 
electroweak corrections are included in the parameter 5'ew = 1.0232 [29]. 

The pole parametrization of the vector form factor is f+{t) = /+(0)[My/ {My— 
t)]. Expanding this expression to second order gives A'| = 2A'^, where and 
are the linear and quadratic slopes. 



hit) = /+(0) 



1 + A' ^ + 



\1 



mi+ 2 

IT' 



We evaluate the phase space integral from the value of My from KLOE, My = 
870.0 ± 9.2 MeV [22], and get = 0.10320 ± 0.00020. The pole fit result is less 
affected by the strong correlation between the linear and quadratic slopes and 
provides better consistency among the values of Ik from different experiments 
(KLOE [22], KTeV [19], ISTRA+ [20,30], and NA48 [21]) than is obtained 
using the results for A'^ and A'|. 

We obtain: 

ff'^'iO) X Vus = 0.2150 ± 0.0014 (19) 

Using f^°'^ (0) = 0.961 ± 0.008 from Ref. 31 (this value is in agreement with 
a recent lattice calculation [32]), we get 

= 0.2238 ±0.0024. (20) 

To perform a test of first-row CKM unitarity, we define: 

^ ^ud ^us ^ub 
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Using Vud = 0.97377±0. 00027 from Ref. 33 and including the small contribution 
of Vub [13], we obtain 

A = (1.7± 1.2) X 10"^ (21) 
which is compatible (1.4 a) with zero. 

A new measurement of R has recently been made at KLOE, with a precision 
improved by more than a factor of two with respect to Eq. 7. Combining the 
new and old KLOE measurements, we obtain R= 2.2549 ±0.0054. The results 
presented in Eqs. 10, 12, 13, 19, 20, and 21 depend only slightly on the value 
used for R. The results obtained using the updated value of R are listed below: 

BR{Ks TT+TT") = (69.196 ± 0.051) x 10"^ 
BR{Ks ttV) = (30.687 ± 0.051) x 10"^ 
BR(Ks TT'e+iy) = (3.528 ± 0.062) x 10"^ (22) 
BR{Ks n+e-p) = (3.517 ± 0.058) x 10"^ 
BR(Ks neu) = (7.046 ± 0.091) x 10"^ 
Prom the last of these, the following quantities follow: 



Rea;+ = (-0.5 ±3.6) x 10"=^ 

ff'^'iO) X Vus = 0.2153 ± 0.0014 
Vus = 0.2240 ± 0.0024 



(23) 



A = (1.6 ± 1.2) X 10 



-3 



The differences between these quantities when calculated using the old and 
new values of R are well within the stated uncertainties. 
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